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Abstract 

Let /i be a probability measure (or corresponding random variable) such that all moments 
exist. Knowledge of the moments is not sufficient to determine infinite divisibility of the 
measure; we show also that infinitely divisible, and in particular lognormal, distributions lose 
infinitely divisibilty when censored in certain ways even if all moments are arbitrarily close to 
those of the uncensored distribution. The moments of a composition of k copies of ^ are ex- 
pressed as combinatorial compositions of the We express the moments of the compositions 
in the context of occupancy problems, arranging n balls in k cells; the classical convolution 
is described by Maxwell-Boltzmann statistics and is multinomial. For certain non-infinitely 
divisible measures with moments increasing fast enough the indexing of a k-cell combinatorial 
composition is extended to indexing by non-negative real t and we construct classical convo- 
lution measure semigroups from amongst the t-indexed classes. We prove also that when a 
random variable with infinitely divisible distribution is embedded in a Levy process (It) then 
the t-indexed Maxwell-Boltzmann is the law of 1^. In order to get moment-based multinomial 
compositions indexed by a continuum we use random measures and random distributions rather 
than random variables. An alternative approach to embeddability of a non-infinitely divisible 
^ is by considering non-classical convolution measure semigroups; for example embedding /x 
in a Boolean convolution measure semigroup and retaining the multinomial character of the 
moments. Embedding /i in an Urbanik generalised convolution measure semigroup one loses 
the multinomial character of the moments. 
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1 Semigroups and non-infinitely divisible distribu- 
tions 



1.1 Moments of censored lognormal distributions 

The (a, a)- /og'norma/ distribution X^^o- has probabihty density 

^jT^-^ ^'^"2^2'"^ }; support R+ and is known to be infinitely divisible with moments 

= e"°+5'^^'^^. A third parameter could describe shift of the support along R. There is 
no moment generating function since the moment expansion of the characteristic function 
does not converge at all. The characteristic function is known ([SOI) and the density 
function for its Levy measure can be expressed as e~^^U{^)d^ ([^); but U cannot 
be explicitly determined. 

For a real random variable X the spectrum cr[X) of X is taken to be the set of all 
points of increase of its distribution function. It is well-known that the spectrum of an 
infinitely divisible distribution fi is the closure of vector sum of n copies of the spectrum 
of /i", for every n G N. 

For any lognormal distribution one can produce a discrete distribution with the same 
moments which is not infinitely divisible. Indeed, it is shown in ^3] that the discrete 

12 2 2 

distribution with weights a""'e~2" °" at the points ae""" ,n G Z has the same moments 
as a lognormal distribution; the parameter a is dependent on a. This distribution is not 
infinitely divisible; indeed, the spectrum of fi cannot be the vector sum of n copies of 
for any n > 1. 

We shall generalise the following: 

Proposition 1.1. If X is a Poisson random variable then P[X G (a, 6)] > implies that 
Y*[X G {na.nh)] > {and also that 

a{X) n (a, 6) ^ implies a{X) n (na, nh) 7^ 0} for alln>l. 

Proof. Obvious since P[X = ra] = e~^X^/n\ for some A > 0. □ 

Definition 1.2. The left-truncated lognormal is constructed by truncating the mass over 
an interval {0,b),b > and transferring the removed mass to the origin. The right 
truncated lognormal is constructed by truncating the interval [c, +00) for some c > and 
redistributing the removed mass over [0, c). The gap-censored lognormal is constructed by 
cutting off the mass over the interval {a,b),0 < a < b, and transferring that mass to the 
origin. 

Proposition 1.3. Right-truncated lognormal distributions are not infinitely divisible. 
Proof. They cannot be infinitely divisible as they have finite support ( |T2j §8.4). □ 
Definition 1.4. For a measure ^ on R+ for which all the moments are finite we call 

2 

C = limsupn — — — the critical ratio for the moment sequence of a. 

In Proposition 14.151 the critical ratio is shown to provide a criterion for indeterminacy 
of the moment sequence (/x^). We show that censored lognormals can produce non- 
infinitely divisible distributions with the same critical ratio as a lognormal. 
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Proposition 1.5. For large n, by choosing small enough b > the moment m„ for the 
left-truncated lognormal distribution can be arbitrarily close to the untruncated nin and 
the critical value C for the truncated moments will be arbitrarily close to C . The same 
is true for gap-censored lognormal distributions. 

Proof. For the purposes of the Proposition there is no loss of generahty in assuming a = 
0, cr = 1. Denote e~^^^ du by \I>(x). Using one sees that ^ = \l>(log6 - n)/'^{\ogh). 
Thus ^ is, for small enough 6, arbitrarily close to 1 and converges uniformly to 1 as 

~ 2 2 

n +00, so lim^^o lim sup„ ^ — ^ = limsup„ — — . The Proposition follows. □ 

Definition 1.6. We call a Poisson distribution with a random intensity mixed Poisson 
and we denote the mixed Poisson with intensity X by X . 

It is well-known that X is infinitely divisible if X is infinitely divisible. 

An infinitely divisible random variable X on Nq has spectrum in N and probability 
generating function '^PnZ'^ = e~^^^~'^^^^^ where p„ = P[X = n] and \Q{z) is also a 
probability generating function, written as Yl'o' ^kZ^ with > for all k. 

Definition 1.7. Katti's test \14i on a discrete distribution on No comprises solving the 
equations {n + l)p„+i = ^i=oPj-krk, j ^ Nq for the given pni a negative implies that 
the distribution is not infinitely divisible. Conversely, the distribution is infinitely divisible 
if all the are non-negative. 

Proposition 1.8. The left-truncated lognormal distribution is not always infinitely di- 
visible; in particular it is not infinitely divisible when a = 0, 
a = 1,6 = -1.4. 

Proof. We do this by showing that for the truncated distribution X the discrete dis- 
tribution X is not infinitely divisible. Consider a Poisson distribution with a trun- 
cated lognormal intensity, so the probability generating function is determined by the 

Pk = 7T^s= ^e^^'^""'^ dx, and the truncated mass f '^''^ e~^^^'^^ is added to vn. 
Applying Katti's test we obtain negative values of for sufficienly large k. 

It was more economical for computation to use N copies of the truncated log-normal 

setting pfc = /i^fe Independent computations on Maple and Mathe- 

matica, for an accuracy of 16 digits and = 10, produced negative rg and rg. □ 

Katti's idea was extended, in [201 , to develop a general criterion for infinite divisibihty 
of a distribution on the half-line. We use this to prove: 

Theorem 1.9. The general left-truncated Xa,o- is not infinitely divisible. 

Proof. A probability measure /i with support in R_|_ is not infinitely divisible if for some 
a; > 

/ ydM= fi{[0,x-y])ydu{y) +-ffx{[0,x]) 

J[0,x] J{0,x] 

for all 7 > and all measures u such that J^^ ^^-^ mm(l, x)di'{x) < +oo 

( pOj Theorem 51.1). Choosing x < b then the left hand side of the equation is zero 
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whilst the right hand side is non-zero. Indeed, /i[0,x]) = yu([0,a; — y\) is the mass which 
had been transferred from the truncated section to the origin, so 7/i([0, x]) > and also 

/(o,x] yd^iv) > 0- ° 

Remark 1.10. Another method of proving non-infinite divisibility of distributions on the 
half-line is by determining the thickness of the tail. For example, truncate the left half of 
the Gaussian distribution on R; by adjusting the mass by moving the truncated mass, say 
to the origin, or by redistributing the left- hand- side mass over the right- hand- side. The 
tail remains asymptotically the same as before adjustment and so will be too thin for the 
distribution to be an infintely divisible (see Theorem 26.1 of ]E3i)- 

Definition 1.11. We define a compound Poisson process to be of the form X{t) = 
Yj where the Yj are i.i.d. non-negative random variables representing the length of 
jumps taking place at the epochs of a homogeneous Poisson process N, and N{t) represents 
the number of epochs which have taken place by time t. In the literature the distribution 
has also been referred to as being completely random or as a mixed or generalised Poisson 
distribution. 

The characteristic function is known to be of the form 
^xit) = exp{J (e*^^ — l)(iz/(x)} where its Levy measure u satisfies 

J min{l,x)dv{x) < oo. The usual correction term in the expression is unnecessary since 
the integral converges without it. 

It is well-known that a distribution on R+ with an atom at is infinitely divisible if 
and only if it is compound Poisson. 

We shall usually interpret X as the distribution for X{1). 

Corollary 1.12. If the jumps of a compound Poisson distribution X are Poisson-distributed 

random variables its Levy measure is discrete, and 

P[X E (a, b)] > implies that P[X E (na, nb)] > for all n > I 

Proof. This is a direct consequence of Proposition 11.11 □ 

Definition 1.13. For X{t) be a compound Poisson process with Levy measure v and 
e > we define X^ to be the restriction of X to jumps of length exceeding e. 

Proposition 1.14. ^xX^) has the form exp{t {e^^^ — Vjdu^x)}. 

Proof. Partition (e, +oo) as IJj=i where the Bj are disjoint Borel sets. Define X{Bk, t) 
to be ^^go.(x) XBki^i''') ~ ) 'vi^-! the sum of the jumps of size in B^. up to time 

t. For fixed t the X{Bk,t) are mutually independent and have characteristic functionals 
of the form exp{t J^^ (e*^^ — l)di'{x)}. The Proposition follows. □ 

Lemma 1.15. P[\X — X^\] > i] > converges to as e ^ 0-\-. 

Proof. Since as e 0+, by a standard argument |X — X^l — in probability. 

□ 

Definition 1.16. A random process Z on (a product of complete separable metric spaces) 
X is called a marked point process if, for any Borel set B iny (the space of 'marks'), 
the number of points in the marginal process Z{B) is finite. 
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We shall show that is a marked point process, the marks given by independent 
heights of the jumps. We need the following: 

If e' < e then cr(Xe) C (y{X^i) and the mass at any point of (j{X^i) is larger than 
the mass at a{X^. By Lemma 11.151 for small enough e > one has [a{X^) fl (a, 6)] C 
[u{X) n (a, h)] and the mass at any point of (y{X^i) remains finite. 

Proposition 1.17. Suppose X is compound Poisson with Levy measure v and choose 
e > 0. Let be a Poisson random variable with intensity u{[e, +oo)). Let denote the 
restriction off to [e, +cxo). Marks are given by mutually independent heights {hj),{j = 
1, • • • , N^}, the jumps with distribution ^^^'^^■^ . Then X^ is equal in law to X]j=i ^j- 

Lemma 1.18. Let {a,b) be a bounded interval in (0, +oo). // 

P[X^ G {a,b)] > then for any integer n > 1 there is an arbitrarily small rj > such 
that also P[X, G {n{a - r]),n{b + r])] > 0. 

Proof. There is a finite set Ji of integers j such that P[Ylj(zj hj G A] > 0. Fixing n, since 
the underlying process is Poisson, for every hj : j G Ji there are also n different jumps 
with heights hk arbitrarily close to hj] we denote these heights by hj,j G J/3, /? = 2, 3, .., n. 
Thus there will be an arbitrarily small 77 > such that P[J2j&jiu uj„ ^ (^'^ ~ ^^^5 + 
nr])] > 0. □ 

Theorem 1.19. Suppose < a < b and X compound Poisson. Then P[X G (a, 6)] > 
implies that P[X G {na,nb)] > for any integer n>l. 

Proof. Choose rj > 0. If P[X G (a, b)] > then P[X^ G [a — rj, b+rj)] > for small enough 
e. By Lemma II.8I one has P[X^ G {n{a — ri),n{b + rj)] > 0. Since X — X^ and X^ are 
mutually independent, using Corollarv 11.121 and Lemma Fl.lHI P[X G {n{a — rj) , n{b + rj))] 
> P[Xe G {n{a - ri),n{b + r])] x P[\Xe - X\ < 77]) > 0. Letting r] 0+ the Theorem 
follows. □ 

Corollary 1.20. Gap-censored compound Poisson distributions are not 
infinitely divisible. 

Proof. Suppose there is a single gap between a and b. If 2a < 6 then the interval (a, 2a) 
supports mass, a contradiction. If 2a > 6 then the interval (6, 2a) supports mass, if 2a = 6 
then the interval (a/2, a) supports mass and again there are contradictions. □ 

1.2 Combinatorial compositions and their convolution semigroups 

The convolution product of copies of a measure gives rise to a discrete convolution semi- 
group over Nq. We show that we cannot embed the measure in a classical convolution 
semigroup over R+. 

The n'th moment of a convolution product of copies of a probability measure is repre- 
sented by E[(^^^-^ ^i)""] where the Xj are independent identically distributed real random 
variables with distribution /x. We admit non-classical interpretations of (^^L^^ Xi)^. 

Definition 1.21. For probability measures such that all moments exist, fi o u defines 
the class of measures having n'th moment 'Yyj=oCffJ'j'^n--j- We call this the binomial 
composition of /i and v. It need not define a unique measure. We extend this idea to 



5 



various combinatorial composition classes of k copies of fi (not necessarily determining a 
unique measure), written say as . We denote the classical convolution of k copies of 
fi by fi''. We express 

{xi + ..+Xk)"' o,sJ2ni+ +nfc=n '^i' ••' whcrc {ti] Tii, . . , Uk) dcnotcs the number 
of ways of arranging n objects, say balls or particles, into k cells such that there are uj 
balls in the j'th cell. The combinatorial composition is called multinomial if 

n\ 

(n;ni, ..,ng) = — 

nil ■ ■ ■ Ugl 

Definition 1.22. A sequence of positive real numbers is called Stieltjes if it is the moment 
sequence of some /i with support in R_|_, and Stieltjes determinate if this measure is 
uniquely defined by the sequence. 

Theorem 1.23. Suppose there is a combinatorial composition fi°'',k G N of copies of fi 
which can be extended to be indexed by t & R+ in such a way that for each t the sequence 
//°* is a Stieltjes moment sequence, continuous at t for each n. Then there is a (classical) 
convolution measure semigroup A{t) such that A{t) has the same moments as though 
A(l) 7^ fi when /i is not infinitely divisible. 

Proof. There is no loss of generality in assuming t G (0, 1). For a fixed t, since is, 
by hypothesis, a Stieltjes moment sequence, there exists a measure, which we denote by 
A(t), with moment sequence If the moment-sequence is determinate then X(t) = /i. 
The theorem follows from the following lemmas and the fact that if fi is not infinitely 
divisible it is not stable and so not embeddable in the semi-group. □ 

Lemma 1.24. Fix t G R+ and denote the integral part of nt by [nt\. Given a sequence of 
measures (A(-)f"'*l) there is a subsequence {n^jl') of N such that the A(-^)["**l converges 

weakly to a measure, say A{t), with the same moments as fi°^ as n oo. 

Proof. This is standard, using Prokhorov's Theorem and the Chebyshev inequality. □ 

Lemma 1.25. The A{t),t G R+ form a semigroup. 

Proof. In order to introduce sufficiently many subsequences of sequences of integers for 
all s, t G R+ we start with s and t rational. Then 

hm A(— * lim A(— * {Sa/n,))' 

k Uk k Uk 

converges to A(s) * A(t) for some subsequence of the sequences [n^^^] and (n^*^). The 
factor where p can be 1 or 0, is there to adjust the expression whenever nt ^ [nt]. 

By hypothesis /i^* is continuous in a neighbourhood of t so A(t) * A(s) = A(s + t) also for 
s,tGR+. □ 

Remark 1.26. If p is not infinitely divisible and the sequence (/i„) is determinate then 
for any combinatorial composition the cannot be a Stieltjes moment sequence for all 
t G R^. It was shown in that the Maxwell- Boltzmann (/x^*), as defined the next section, 
form Stieltjes moment sequences if all the — — — are all smaller than about . 
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2 Combinatorial convolution semigroups 



2.1 Maxwell-Boltzmann statistics 

We investigate combinatorial convolutions interpreting the n'th moment of combinatorial 
convolutions fi°^ as weighted combinations of the arrangements of n balls into k cells so 
that g = 0, 1, ■ ■ ■ of these are empty. To determine the weights one needs some kind of 
statistics; there is no non-probabilistic formulation. Differences arise if one distinguishes, 
or not, between the various particles labelled xi,X2,--- and/or distinguishes, or not, 
between the cells (say by selecting their order). 

Consider a mechanical system of n balls and k cells such that each ball is assigned 
to a cell. A state of the entire system is described by a random distribution of the 
particles in the cells. We can distinguish between different phase-spaces, for example, 
with Maxwell-Boltzmann, Bose-Einstein or Fermi-Dirac statistics (see §3). To apply the 
above to moments of measures we identify a cell with a measure fi, and q balls in fi is 
represented as fjjq . 

There fc" possible ways of distributing n distinguishable balls in k distinguishable 
cells. For the Maxwell-Boltzmann model, originally proposed to explain the distribu- 
tion of sub-atomic particles into different energy states, both the balls and the cells are 
distinguishable and there are ^ii-'-nki s^^^^sj given cells numbered 1, ■ ■ ■ ,n and allowing 
multiple occupation and empty cells this is the number of arrangements such that nj 
balls are in the j'th cell. Each of these states has the same probability of occurring i.e.. 

Definition 2.1. 03^ = ^^^Zq (-l)-''Cj=(fc - j)", zero forn>k>0. See II We 
shall call these Boltzmann numbers. 

There are 03^ distributions leaving no cells empty and the number of distributions 
leaving q cells empty is So also C^_^^" is the number of distributions with only 

i cells occupied. 

Other characterisations of !B. 

= Eni=i ■ ■ ■ Y2.,=i X(ni+-+n^=n) „,^".'„^.p whcrc X siguifies the characteristic function 
of a set (cf. [iHj §61. 

= j!©" where denotes the Stirling subset number (i.e., the number of ways of 
partitioning n objects in j non-empty sets, identical to the Stirling number of the second 
kind) . 

03" = A-^xj^^Q, where A-^ represents the j'th power of the finite difference operator (cf. 
[nj §2.3). 

Definition 2.2. Given fi, and A; G N, we define the Maxwell-Boltzmann fi'^ as 

k n n ! 

J2^jJ2"'Yl X(ni+...+n,=n) , ^'^^ 

j=l ni=l "j=l i=l 

For t not an integer Cj denotes the usual extension of for the expansion of (1 + x)*. 



7 



Since C = when j > k we can define the Maxwell-Boltzmann /i° as 



Proposition 2.3. T/ie Maxwell-Boltzmann arrangement is multinomial. 
Proof. For A; G N the Maxwell-Boltzmann-/i°^ is indeed 

^-^ nil ■ ■ ■ Uk'. 

711 + .. +nfc=n *=1 

as the expression in the definition is only a re-arrangement of the terms in this sum. □ 

The moments are well-defined, involving finite sums of finite products of moments 
independent of t. The above definition can be confusing for computations. Seeing that 
the 6^ are well tabulated it is more convenient to compute as the polynomial 

n n j 

We elaborate on this in §4. Note that 23" = 1 and fB^ = n\ for all n and 03^ = for 
k > n. 

Infinitely divisible distributions. 

Since we are considering positive random variables, in order to simplify notation we 
shall occasionally use the (possibly formal) moment generating function instead of the 
characteristic function. Assuming all moments of fi are finite, the cumulant generating 
function 6 i— > "^{0), i.e., the log of the moment generating function, can be expanded as 
the infinite series 

Ki9-\-^{9y + - ■ ■ where /tj is the i'th cumulant, viz. — . One compares this expansion 
with the infinite series obtained by expanding the characteristic function into powers of 
6; by using de Faa's formula for derivatives of a functional one obtains expressions for 
the fin in as polynomials in ■ ■ ■/€„) (see [Ej §2.4). 

Definition 2.4. Let Y be a real random variable with infinitely divisible distribution /i. 
We may assume that the characteristic function ofY is of the form exp{i9a + /o°°(e*^^ — 
l)di'{z)}. For t G R+ let Yt denote the random variable with characteristic function 
exp{t[iea + (e^^^ - l)dp{z)]}. 

The characteristic function of an infinitely divisible jj, is of the form e"^*-^-*, so simplifying 
de Faa's expansion. One can show, using Taylor and McLaurin expansions around 6* = 0, 
that 

(e*(^))=n! (^(0)("i)---^(0)('="^))) 

S=(^ mH \-kn^.=n 

giving an expression of moments in terms of cumulants. The cumulant generating func- 
tional for Yt is thus 



d9' 
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and the moments thus obtained are polynomials in t with coefficients polynomials in the 

Theorem 2.5. For infinitely divisibleY with distribution fi the classical Maxwell-Boltzmann- 
fi°^ is the law of the Levy process Yt. 

Proof. The replacement of the symbol A; by t in the Maxwell-Boltzmann and the 
insertion of the multiplyer t before the cumulant generating function of Y to get the 
cumulant generating function of Yt are essentially trivial. Both and E[(yt)"] are 
polynomials in t vanishing at integer values of t, and as such are identical. □ 

Illustration. 

The simplest example is that of a Poisson random variable Y with distribution and 
mean A. All the cumulants equal A. Using the difference-operator formulation of 03 one 
can show that /i„ = Yl^j=i ®i ff ^^"^ hence 

n{YtT] = E;=i As /ii = A, /i2 = A + X\ = A + 3A2 + A3, • • • 

so E[(F()] = tA, E[(Ft)^] = At + (At)^, ■ ■ ■ For t of the form \ one has, for example, 
E[Fi] = i^i, E[(r 1 )2] = - \{^^^)\ Y.[{Y.f] = i/ig - |/i2/ii + ■ ■ ■ • E[y.] = 

Ifii, E[(Fi)2] = |/i2 - E[(Fi)3] = |/i3 - |/i2/ii + ■ ■ ■ , which agree with 

the corresponding combinatorial expressions for the Maxwell-Boltzmann-/in'' and these 
could be used to construct the expectations for Yt with t G Q 

2.2 A continuum of cells 

Since n is a fortiori not stable it is not possible to embed fi in any of these convolution 
semigroups or to interpret this as as an arrangement of objects in a continuum of cells (say 
in a lattice with spacing decreasing to zero). For a chain ... ^ of random variables 

E[(Ei=i ^iT] is well-defined. When the Xi are i.i.d. we have linked it to arrangements of 
n elements in k cells. A limiting continuum of cells would be in the form E[(^^gjQ X^)"] 
but this can have a meaning only if at most a finite number of Xg are non-zero. However 
by taking X to be a random measure or random distribution the continuum of cells can 
have a meaning. 

Definition 2.6. A random measure is taken to be a mapping associating a real random 
variable Xa to each element A of a family of Borel sets in R+j or associating a real 
random variable X{(f)) to test functions cj) which are Borel on R_|_. A random (Schwartz) 
distribution associates random variables to test functions 9 which are infinitely differen- 
tiable with compact support. 

Following [Hj, a random measure X is said to be weakly decomposable if, for every se- 
quence of disjoint Borel subsets ^1,^2,...., the X^- are mutually independent and ^X^- 
converges almost everywhere to Xy^„. The distribution of a weakly decomposable X is 
interpreted as the continuous product of the distributions of its components It follows 
from jH] Theorem 2.1 that for a weakly decomposable X each Xa is infinitely divisible 
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and ErG(s,t)^M = ^(s,t)- 



If X^tj is zero for each t then the characteristic functional is of the classical Levy- 
Khinchin form. 

A function / is called decomposable if for every finite partition Ai, - ■ ■ ,An of [0, 1] 
into measurable sets / may be written as /i + ■ ■ ■ + /„, where each fj is measurable with 
respect to the field of measurable sets generated by the Aj. The conjecture in |H] that, 
for the (T-field of Borel sets in [0, 1], there are enough decomposable functions to generate 
all measurable functions was verified in [221 ^y showing that t ^io,t) is a Levy process. 
The characteristic functional factors into a product of independent Gaussian parts of 
various dimensions and a Poissonian part (see [H] Theorem 4.1). 

What kind of processes, denoted say by Y{t), can be constructed from random mea- 
sures and distributions? One such process is t \—>- Y{t) = Xp^t) as above. For random 
distributions, if t t-^ M{t, 6) is integrable with respect to t and 6 e^o M{s,es)ds jg pQgj_ 
tive definite (to ensure that it is the Fourier transform of a Schwartz distribution) then 
eIo^Hs,8s)ds ^YiQ characteristic function of X]se[o ^ (Gaussian) white noise 

corresponds to the characteristic functional with M(t,9) = —\(P'- Also a Poisson-type 
continuous product can be constructed based on a Fermi-Dirac system (cf.|^). 

These above are Focfc, also called linearisable, continuous products. The more inter- 
esting non-linearisable continuous products have been investigated in |2S| using partitions 
other than the Borel cx-fields. 

Moments. To get moments of the continuous product we could compute E[(Xi)"|Xs] 
or E[{Xty'Xs{(f))] where say X(0) = EfcLo^fc^'' = ^''^ ■ Matter case for 

t G [0, s] we may write 

E[(Xt)"e*^^i] = E[(Xt)"e*^^=]E[e^^(^=~^*)]. When E[e^^^»] is of the form e'^^'^ one can 
calulate the moments by differentiation of the 

characteristic functional. However we already know by Theorem 12.51 that when Yt con- 
structed from X is a Levy process the moments E[(y()"|Yi] are multinomial. 

3 Non-classical convolution measure semigroups 
3.1 The generalized convolution 

A generalised convolution introduced by K.Urbanik, is a binary operation, indicated here 
by ■, on non- negative probability distributions, associative commutative and distributive 
for convex sums. The translation mapping by x G R+ is denoted by Tj,. It is assumed 
also that there exists a sequence of real numbers c„ such that Tc„6'{^ converges weakly to 
a measure cTk ^ 6q, called the characteristic measure. Several examples are given in [TTj . 

The convolution is called regular if it admits a (generalised) characteristic function 
such that there is a bijective correspondence between probability measures and charac- 
teristic functions. The translation mapping by x G R+ is denoted by T^. It is postulated 
that there exists a sequence of real numbers Cn such that Tc„6'{^ converges weakly to a 
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characteristic measure a^. ^ Sq. The characteristic function of is °° exp~^"'^dai^{x) 
where k is called the characteristic index. The semigroups defined by the generalised 
convolutions correspond to those for ■ -Levy processes as described in |25j . 

For the classical convolution = 5o and k = 1 and the characteristic function is the 
Laplace transform of the measure. 

Proposition 3.1. No generalised convolution, other than the classical convolution, is 
multinomial. 

Proof. For K = 1 and k > n, /i„^ = YJ'r=i (-1)"^''^^ C'n~r^(/^n), as in |2HI. This is 
compatible with the multinomial combination. For k 7^ 1 the expansion is incompatible 
with the multinomial expansion. □ 

3.2 Non- commutative systems 

A non-commutative probability theory can be constructed as a non-commutative algebra 
A of observables with a normalised positive linear state, say 0; when A is commutative to 
X & A there corresponds a measure such that fJ'x{4>) = (pix). Bose-Einstein and Fermi- 
Dirac statistics relate to non-commutative quantum systems. Combinatorial convolution 
semigroups of measures derived from quantum 'random variables' relate to consideration 
of operators on Fock spaces. Eliminating the (Hamiltinian-based) dynamics and fixing 
the number of particles, so ignoring creation and annihilation, one obtains a commutative 
subalgebra of A generated by the number-operators. For Bose-Einstein and Fermi-Dirac 
systems the cells correspond to energy levels, 'balls' are indistinguishable and only dis- 
tinguishable arrangements are considered so (n; ni, .., n^) = 1. For Bosons there will be 
Cn^^~^ states. For Fermi-Dirac statistics, by Pauli's exclusion principle, a cell will be 
either empty or contain one ball; there are less balls than cells so there will be states; 
in the computation of (X]^=i^«)" o'^^ postulates XiXj = —XjXi for i 7^ j. Convolu- 
tions of Bose-Einstein and Fermi-Dirac random variables have been constructed in [23 
by consideration of graded vector spaces and bi-algebras. The Fermi-Dirac convolution 
can represent a one-dimensional spin | Heisenberg ferromagnet (the ferro-magnetism as 
proposed by Heisenberg) sitting in a Fock space, whereas the Bose-Einstein model has 
integer spin. Quasi-particles, called magnons, carry the magnetic domain field in this 
model. 

Remark 3.2. Seeing that a continuum ferromagnet does not exist (see |^ and \24]l ) there 
is no chance that the Fermi-Dirac convolution product fi°" can be extended to a as we 
do for Maxwell-Boltzmann statistics. 

3.3 The Boolean convolution 

The Boolean convolution product, related to a non- commutative probability theory, is 
denoted by l±) has been dealt with in [22] • A measure fi can be identified, by its moment 
sequence, with a (normalised positive linear) states on C < X > and /il±)i/ is identified 
with states on C < Xi, >, where C < Xi, X2, X„ > denotes the ring, with complex 
coefficients, of polynomials in noncommuting indeterminates Xi, X2.., X„. The restriction 
in [22] that /x has compact support is unnecessary and the results are valid under the 
condition that all moments exist (but then a moment sequence need not be determining). 
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In jl] the Boolean convolution semigroup R+ 3 t fi^* was shown to exist. The 
semigroup is not the same as the classical convolution semigroup; the characteristic func- 
tion of is a 'self-energy' based on the non-linear Cauchy transform. It can be expressed 
as a -f Jj^^ ^j^dr^z) for a e R and some measure r on R, this form being conducive to 
defining infinite divisibility of arbitrary distributions. Indeed it is shown in j22] that all 
/i, as above, are 'infinitely divisible'. 

The Boolean convolution is related to a lattice of non- crossing partitions; these are 
used also in quantum-type convolutions (cf. §2), 

Proposition 3.3. The Boolean convolution is not a bilinear relation; indeed it is highly 
non-linear and it is not one of Urbanik's generalised convolutions. The Boolean convolu- 
tion is binomial but not multinomial. However the moments of the fi^^ satisfy multinomial 
decompositions and /i^^ = /i^ for n, A; G N. 

Proof. It not distributive for convex sums; for example, |(5a W \Sb) W (5c is a convex sum 
of ^a+b+c±^y/ {a-b)^+c'2 ^^^^ Example 3). For n > 2 the measure /i'*'* does not have the 
same n'th moments as any element of For example, for random variables X, Y, Z one 
has 

(/i W u){XYX) = i2{X)iy{Y)fi{X) so 

{fi l+l z/)3 =1713 + 27712^1 + 2n2mi + (mi)^ni + (ni)^mi + ns 

as opposed to the classical (/i * z/)3 = m-^ + 3m2?2i + 3min2 + n^. The remainder of the 
Proposition follows as the fj!^^, A; G N, are mutually commutative. □ 



4 Log-convex sequences, Stieltjes moment sequences 
and total positivity 



and denote it by On- We shall call the sequence 6-log-convex if 9n < for all n. We call 
a 6-log-convex sequence strictly log-convex when 6 < 1. We similarly index log-concave 
sequences. 

We may assume, without loss of generality, that Oq = 1. The sequences of moment 
ratios were used by T.J. Stieltjes in ^| Chapter II. 

Remark 4.2. a). It is shown in Theorem 51.3 of 120^ . using Katti's test that a 
probability measure on N is infinitely divisible if (pn) is log-convex, 
b). Log-convex and log-concave sequences are known to be unimodal in the sense that they 
are initially decreasing (increasing) with n until one or two nodes uq and no + 1 and then 
increasing (decreasing) with n. 

For (a„), log-convex and log-concave with moment ratios less and greater than 
61,62 respectively, the sequence (anbn) will be log-convex if 6162 < 1. 

The binomial sequences and the sequences {&^)k of Stirling subset numbers are well- 
known to be log-concave. 

Proposition 4.3. Stieltjes moment sequences are log-convex. 



Definition 4.1. Sequences (an) of positive reals such that the 
are called log-convex, or log-concave, respectively. We shall call 





a moment ratio 
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Proof. It is obvious that /i„ is monotone non-decreasing. Since, as shown in ^^^^ is 
monotone non- increasing, ^„ = -tn—IEn+l. < x_ □ 

The following proposition for log-convexity can be adapted also to log-concave se- 
quences. 

Proposition 4.4. Suppose the sequence (an) of positive real numbers has On ^ 6 < 1 for 
all n. Denote log an by fn- 

a) When 6 < 1 one can extend the sequence (fn) to a smooth continuous convex function 
F on the relevant interval o/R+ with F" > log-^J^, 

F{t) > e*''°3v^ + 0{t) and so an > e"''°^^. Here F'{t) diverges as t ^ oo. 

b ) When sup On = 1 the extension of the fn to F is such that 

limj^oo F'{t) = p < oo; F{t) is asymptotically linear so an is asymptotically the order of 

Proof. We construct F by fitting parabolas touching tangents, determined by chords, at 
the points (n, f{n)). The 'acceleration' of F at is 

|(loga„+i + loga„_i - 21oga„) > -logVO. □ 

Remark 4.5. Case (b) cannot occur for Stieltjes moment sequences unless the moment 
sequence is determinate. If the characteristic function is analytic then lim„ sup — — = 

^n— l^n+1 

2 

1; there is large class of such u having lim„ — — = 1 but this does not hold in general 

(see 112]). 

Corollary 4.6. // (a„) is strictly log-convex then for k < n one has 

< I. If (a„) is -log-convex then for k < n one has < Q^^"^-^) , 

Proof. From the concavity of F the gradient of F is non-decreasing as n increases. Con- 
sidering /o, fkfn-k, fn the first part of the Corollary follows. The second inequality holds 

because the sequence (6' 2 a„) will be log-convex. □ 



Proposition 4.7. Fixing n, the sequence {^1)k is log-concave and k \—>- 03^ can be 
extended to a continuous convex function on (0, n) . 

Proof If J < n - 1, then by [IE], > so 

> -rj— > 1. The second part of the Proposition follows from Proposition 



01 □ 

We note that the ^„ — will be close to 1 except for small values of n so the 
sequences 03^ Y[l=i f^ni will almost always be log-convex. 
Definition 4.8. A sequence (aj) is called positive- definite if 

Ylij=o^i^j^i+j — /^^ '^'^'^ '^''^y sequence of non-zero (q), strictly positive- definite 

when the sum is greater than 0, semi-definite if the sum can be zero. 
The Hankel matrix [ao ■ ■ ■ a2n] 'is defined as the (n + 1) x (n + 1) matrix having leading term 
uq and aij = Ui-^j for all i,j < n. Its determinant will he written as A„ when it is clear 
what sequence is involved. The determinant of the Hankel matrix [a^, • • • , ap+2n] will he denoted 
likewise hy A„^p. 
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A necessary and sufficient condition, going back to Stieltjes jSHj, that a sequence 
(a„), a„ > be a Stieltjes moment sequence is that both 

n t— > a2n and n ^ a2n+i are positive-definite functions for all n, and so also if and only if 
both [ao, ■ ■ ■ , a2n] and [ai, ■ ■ ■ a2n+i] are positive definite for all n. 

Definition 4.9. Let M he a square matrix with coefficients say Mij. Let a = {ai, ■ ■ ■ ak), (3 = 
(/?!, ■ ■ ■ Pk) where the a and (3 are made up of non- decreasing sequences of integers. De- 
note by M[a\l3] the minor obtained using rows and columns numbered a and (3 respectively. 
When a and (5 are made up of consecutive integers we call M[a\(3] a Fekete minor. 
The matrix is called positive- definite if ■ CiCjMij > or equivalently if all the principal 
minors are non-negative. 

It is called totally positive if every minor is non-negative. We shall use the term strictly 
totally positive if the above minors are positive, and totally semi-definite if any of the 
minors are null. 

Theorem 4.10. M. Fekete /?/ Satz U. A necessary and sufficient condition for a matrix 
be strictly totally positive is that all Fekete minors are strictly positive. 

Theorem 4.11. Let (a^) be a sequence of positive reals. When (a„) is a Stieltjes moment 
sequence then its associated Hankel matrix is totally positive. If all the Hankel matrices 
[op, ■ ■ ■ , ap+2m] o,re totally positive then (an) is an indeterminate Stieltjes moment se- 
quence. 

Proof. By definition, a totally positive matrix is also positive-definite. For (a„) a Stietjes 
moment sequence both n — > and n — > a2n+i are positive-definite. Any Fekete minor 
with leading term a2m can be identified with a principal minor by considering the sequence 
n a2n with leading term a2m and any Feteke minor with leading term a2m+i with a 
principal minor for the sequence n —>■ a2n+i- Similarly mapping from principal minors for 
either n — >■ a2n or n — >■ a2n+i one can map to all the Fekete minors. Therefore the Hankel 
matrix will be totally positive by Fekete's theorem. □ 

Theorem 4.12. a). A Stieltjes moment sequence is indeterminate if and only if both 
n\-^ an and n t-^ an+i are stricly positive definite J7W (cf. Theorem 2.18 of 121]). 

b) . A Stieltjes moment sequence is indeterminate if and only if both the sequences ( ^^" ^ ) 

and ( ^^^'"^ ) have non- zero limits as n ^ oo 111^ (cf. \l^ Lemma 3). 

c) . Let (fin) be a Stieltjes moment sequence. Consider the 2-dimensional section = 
{{x,y) : {x,y, fj,2, • " )/"2m)} ofthe cone in (2n+l) -dimensional space formed by vectors {fiQ, ■■ ■ ,/i2m)- 
Let ci^m be such that Ai^2m vanishes if one replaces /ii by ci^m- The sequence {ci^m) is non- 
decreasing, bounded by ni and we denote lirrim^oo ci,m by ci. When c\ = fii the sequence is 
determinate. By \1^ Theorem 2, the sequence is indeterminate if and only if /ii > ci and 
(ci,/xi] C D where D = Hm-^m denotes the limit parabolic region as defined there. 

Proposition 4.13. A measure on R+ with moment generating function 
convergent in some neighbourhood of the origin is Stieltjes determinate. However a Stielt- 
jes moment sequence may determine a measure which has no well-defined moment gen- 
erating function. 

Proof. The first part is obvious since the generating function will define the measure 
uniquely. Distribution functions / with determining Stieltjes moment sequences but no 
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well-defined moment generating function are easily found using G.H.Hardy's (1917, well 
before the Carleman) criterion for determinacy viz., for some q > 1 and 6 > one has 



Hardy's criterion can also be adapted for discrete distributions. Other suitable mo- 
ment sequences could be constructed, using Theorem I4.12r c) and manipulating Stieltjes' 
Hankel matrices. 

Remark 4.14. a). The critical index 6 does not provide a necessary condition for a 
log-conves sequence to be Stieltjes or for a Stieltjes moment sequence to be indeterminate. 
Lognormal distributions, as in §i, do not satisfy the 6 -condition if > log-^ > log2 but 
are indeterminate no matter how large the moment ratio. 

b) . By Proposition \4.!^ a moment sequence is always log-convex. The advantage of log- 
convexity is that one needs to compute the strict positivity or semi-definiteness of only 
one of the shifted sequences. 

c) . The authors of [5], studying total-positivity for matrices, not necessarily Hankel, 
get a slightly smaller value for 5. They do not use a general sufficiency condition but 
compute inequalities using selected critical minors and relating the critical index to a rate 
of decrease of minors away from the diagonal. 

d) . Can the critical moment ratio be shown to be j? 

Corollary 4.15. The indeterminacy of a Stieltjes moment sequence depends only the 
asymptotic values of On as n —>■ oo. 

Proof. When one knows that the sequence is a moment sequence, for indeterminacy 
by part (a) or part (b) one need consider only the asymptototic values of the 6n as 
n —>■ oo. □ 

See Proposition 11.51 regarding our use of this Corollary. 

Remark 4.16. The moments of lognormal distributions, described in §i, form totally 
positive Hankel matrices. Indeed the matrices are easily seen to have generalised Vander- 
monde determinants, well-known to be totally positive. 

Corollary 4.17. // (a„) is positive semi-definite every partial sum 
CiCjai^j for i,j<n can be positive for any fixed n. 

Remark 4.18. Sufficient conditions for positive definiteness of a matrix involving the 
off-diagonal M[(z, j)] being small enough are known. For instance, it is not hard to 
prove that if 'Ylii^j j)] < 1 then the matrix is positive definite. This relates 

to the decrease of off-diagonal minors used in In the case of Hankel matrices it is 
simpler to depend on having a fast enough rate of increase of the fin- 

As in PI a sequence of positive elements is a Stieltjes moment sequence if the growth 
rate of the a„ is fast enough. Indeed the growth rate 




u > 0, ^ < a < 1, noting that here the Hardy integrals are Laplace transforms. There is 
no moment generating function since e^^e~'^°' du = oo. □ 
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c^n+i > ((^ + l)!ctn)"^^ suggested there is fast enough for the sequence to be indetermi- 
nate. It is essentially sufficient that the > (I)"- 

The Maxwell-Boltzmann convolution semigroup. 

It will henceforth be sufficient to restrict t to be in (0, 1). The phenomena evident in 
the expansions of /i^, illustrated in §2 for Poisson /i, persist for the Maxwell-Boltzmann 

for a general /x and for this we will be referring below to the the alternative expansion 
in §2 

Lemma 4.19. Being unimodal (see Remark \4-^ the sequence (23^), for fixed n, is at 
first increasing and eventally decreasing. It is known (see fTW), that 

If = (¥)"iTT { 'tl^r^' > ^ (¥)"ITT'- ^^M-^^ 'hat 

< i^^)^ (k+i)2 ■ ^ fixed, is asymptotically ^ so 03)! is asymptotically k^. In 
particular, 2"-^ - 1 < 6^ < 2"^^ so 2" - 2 < < 2" for all n. 

Expanding the series, 

+ C:_,fl^fi^^^) + C*®^(AX-3/^2/ii + ■■■) + ■■■ + Cin\{fi,r- 

The term accompaning ^2 accounts for when there are n — 2 empty cells. The coeffi- 
cients 2" — 

are the probabilities of the various arrangements fXn-jli'j, J = 1, 2, ■ ■ ■ , n—1. 
The is put in the leading position as it will be larger than the other products of 

moments for a reasonable rate of increase. In the term corresponding to more non-empty 
cells the probabilities for the products of three moments are more complicated. We have 
not elaborated on these probabilities, denoted A]., but inequalities between the *B^, anal- 
ogous to those used in Lemma 2 of [Q, can be obtained by using CoroUarv 14 . 61 and Lemma 

Km 

Lemma 4.20. Let (fin) be a strictly log-convex Stieltjes moment sequence. For t G 
(0, 1) each can he written as a polynomial with leading term tfin, subsequent terms 
alternating in sign, because the terms of C* has alternating signs, and decreasing in 
modulus. After any term in the sequence the absolute value of sum of the subsequent 
terms is smaller than the absolute value of the preceeding term. 

Proof. The absolute decrease of the terms is proved using Lemma 2 of with e = 1; 
Corollarv l4.6l gives Hn-kl^k < fJ'n and this can be used to reduce products of k moments to 
products of k-1 moments, and so on. The last part holds generally for alternating series 
with terms decreasing in modulus. □ 

Corollary 4.21. Given a 9-log-convex moment sequence, one has 
tfin > lJ°n > (1 ~ 6')t/i„ for all n. 

Proof. By Lemma I4.2UI we know that < and 

/^°* > t/^n + ^25^32^2 (Efc^J C^/^„_fc/ifc). We use > 2" - 2. From Corollary EH we 
have < 0^('^-'^) so ^^=1 CJ!^n-kIi'k < 26'/x„. The Corollary follows. □ 
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Proposition 4.22. For a 6-log-convex Stieltjes moment sequence {fin) with small enough 
6, or equivalently for a large enough rate of increase in the moments, for any t G R+ the 
Maxwell- Boltzmann is a Stieltjes moment sequence. 

Proof. It is sufficient to show tliat for a large enough rate of increase of the moments one 
has ai^" \t < 5 for all n. We know also that Un+iUn-i > il^n)'^: the latter being the 

leading term in the denominator. After cancelling out the common factor t"^ the leading 
term, in the denominator, will for fast enough increase of the /i„ dwarf all other terms 
so the Jf" L will be small enough to ensure that (yU°*) be an (indeterminate) Stieltjes 
moment sequence. 

Furthermore, by Corollarv l4.21l i^" L < • As in 111 it is sufficient that 77-^ < 

S, from which one can calculate the maximal 6 that ensures that (f^n) is Stieltjes. □ 

Conjecture 4.23. The maximal 9 from the proceeding proposition is slightly less than 
|. Because, as in Corollary \4-^l\ the inequality has been obtained using only the first 
two terms of the expansion of we conjecture that the estimate can be improved by 
consideration of further terms in the expansion of ix°^. 
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